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Abstract 

We present 0{a^) QCD corrections to the fully-differential decay rate of a fe-quark into inclusive 
semileptonic charmless final states. Our calculation provides genuine two-loop QCD corrections, 
beyond the Brodsky-Lepage-Mackenzie (BLM) approximation, to any infra-red safe partonic 
observable that can be probed in & ^ Xuei/ decays. Kinematic cuts that closely match those used 
in experiments can be fully accounted for. To illustrate these points, we compute the non-BLM 
corrections to moments of the hadronic invariant mass and the hadronic energy with cuts on the 
lepton energy and the hadronic invariant mass. Our results remove one of the sources of theoretical 
uncertainty that affect the extraction of the CKM matrix clement | Vut \ from charmless inclusive 
B-decays. 



Studies of CP violation in B-mesons per- 
formed by BELLE and BABAR, firmly es- 
tablished the correctness of the Cabbibo- 
Kobayashi-Maskawa paradigm at the few per- 
cent level. These studies will be continued, 
when the super- i? factory in Japan will come 
on line. A powerful tool to test the CKM pic- 
ture is the unitarity of the CKM-matrix. A 
combined fit to all available data gives \Vub\ — 
3.58(13) X IQ-^ [1] . This number should be com- 
pared with the value \ Vub\ = 3.38(36) x 10~^ ex- 
tracted from exclusive B — hev decays, where 
the hadron h is either a pion or a p-meson, and 
with iKbl = 4.27(38) x 10"^ which is obtained 
from inclusive measurements of i? — XuCv de- 
cays 0. Although exclusive and inclusive re- 
sults are not in serious disagreement, they are 
clearly different and further scrutiny of both 
exclusive and inclusive determinations of \Vub\ 
is certainly warranted. 

The most complicated theoretical issue for 
both exclusive and inclusive methods is the con- 
trol of non-perturbative effects. This is hard to 
do for exclusive decays and important input in 
this case is provided by ab initio lattice QCD 
calculations of exclusive i? tt, p, .. transitions. 
In contrast, in case of inclusive semileptonic 
decays of i?-mesons, non-perturbative difficul- 
ties can be largely circumvented by the applica- 
tion of local operator product expansion (OPE) 
d, II Hi- The OPE allows to compute suffi- 
ciently inclusive observables related to semilep- 



tonic decays of _B-mesons, such as the total 
rate and moments of various kinematic distribu- 
tions, by correcting distributions and rates of 
semileptonic decays of fo-quarks with a limited 
number of universal non-perturbative parame- 
ters. These non-perturbative parameters can 
be determined from fits to semileptonic decays 
of S-mesons to charmed final states B — > XcCD 
[1, 0, i, i, E, mi and then used in the de- 
scription of B — > XuCv transitions, facilitating 
the extraction of the CKM matrix element | Vub \ 
from observables in the latter. 

While this procedure is well-defined theo- 
retically, it was not used in the determination 
of I I right away because B X^ev tran- 
sitions suffer from a much larger B — >■ X(.ev 
background. One can place severe cuts on the 
kinematics of final state particles to suppress 
it; for example, requiring that the hadronic in- 
variant mass is smaller than the mass of the D- 
meson, ttid 1.87 GeV, clearly eliminates the 
charm background. However, it was realized 
early on that such cuts lead to problems with 
the convergence of the operator product ex- 
pansion and infinitely many terms in the OPE 
need to be summed up to obtain reliable re- 
sults. Such a resummation is usually expressed 
through the so-called shape function [l^, 
which parametrizes the residual motion of a 
heavy quark inside a heavy meson. A recent 
discussion of i? — > X^eD decay in the shape- 
function region, that includes next-to-next-to- 
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leading order (NNLO) QCD effects, can be 
found in Ref. jl4|. Unfortunately, current un- 
certainties in the functional form of both lead- 
ing and sub-leading shape-functions are signifi- 
cant and affect a precise determination of | Vub \ ■ 
In parallel to the studies of the shape func- 
tion region, it was suggested that a combina- 
tion of cuts on hadronic and leptonic invari- 
ant masses |15i| allows one to extend the phase- 
space coverage in B ^ X.^eD decays and make 
the impact of the shape functions smaller. Mea- 
surements that use these selection criteria were 
performed by the BELLE collaboration [l6| . 
Further advances in experimental techniques 
allowed to achieve an almost complete phase- 
space coverage in B ^ XuCP decays. Indeed, in 
recent experimental measurements it was pos- 
sible to fully reconstruct the BB kinematics 
from their decay products, thereby allowing to 
extend selection cuts for the b ^ u process into 
the charm-rich regions and yet, successfully re- 
ject the b XcCv background. For example, 
two recent measurements by BELLE and 
BABAR [l^ present partial decay rates and 
a variety of kinematic distributions for b u 
transitions with the cut on the electron energy 
as low as i5; > 1 GeV. These cuts are inclusive 
enough so that the local OPE expansion can 
be used with confidence to describe B Xuev 
decays. 

We summarize now the status of the theo- 
retical description of _B — > X^eD decays, under 
the assumption that the local OPE is applica- 
ble. The OPE expansion in the inverse 6-quark 
mass rrib is well-established for moments of the 
hadronic invariant mass and the hadronic en- 
ergy d, 0, 0]. The leading order term in the 
OPE expansion is given by the partonic b ^ u 
transition. The total decay rate for b ^ u 
is known in perturbative QCD through ©(a^) 
[igj and a large number of kinematic distribu- 
tions and their moments are known through 
0{as) 0,[2i|,|2i,|2i]. Also, the so-called BLM 
©(^oCKg) corrections that can be derived 
by considering the contribution of a massless 
qq pair to the b u transition, are known 
for the decay rate and main kinematic distri- 
butions 25|, |26j, |27|. The only kinematic dis- 
tribution in 5 — )■ M decays that is known be- 
yond the BLM approximation is the electron- 
neutrino invariant mass distribution, computed 
in Ref. [1^. While the BLM-approximation is 
known to account for a significant fraction of 
the complete 0{a1) correction, the precision 



of current and, especially, forthcoming mea- 
surements of I Vub I , the relatively large value of 
as{mb) and a large variety of kinematic cuts 
employed in experimental analyses make it very 
desirable to compute NNLO QCD corrections 
to the fully-differential 6 — w decay rate be- 
yond the BLM approximation. The goal of this 
paper is to provide such a computation. 

The calculation of NNLO QCD corrections 
to the b — >■ ueD decay requires three ingredients: 
i) two-loop amplitudes for the 6 — > uev transi- 
tion; ii) one-loop amplitudes for b — > ugeD; in) 
tree amplitudes for b — > uggeu and b uqqev. 
The two-loop amplitudes were computed by 

several authors in recent years m (sain, nil. 

The one-loop amplitudes for b — > ugev can 
be extracted from the computation reported 
in Ref. [3^. Finally, the tree amplitudes for 
b — >■ uggeu and b — >■ uqqev are straightforward 
to calculate and compact results can be ob- 
tained by using the spinor-helicity formalism. 
These amplitudes can be found in Ref. 34]. 

The well-known challenge for fully- 
differential NNLO QCD computations is to 
put these different contributions together in a 
consistent way. This is not easy to do since 
individual contributions exhibit infra-red and 
colhnear divergences and correspond to pro- 
cesses with different final-state multiplicities. 
For the computation reported in this paper, 
we use a method proposed in Refs. 35|, |36| 
(see also [37|) which combines the idea of 
sector decomposition fssl Isol l40 l| with the 
phase-space partitioning 41| in such a way 
that singularities are extracted from matrix 
elements in a process-independent way. This 
framework leads to a parton level integrator 
which can be used to compute an arbitrary 
number of kinematic distributions in a sin- 
gle run of the program. We have recently 
given a detailed description of the relevant 
computational techniques in a paper (33 | 
that describes a calculation of NNLO QCD 
corrections to a related process t — > be'^v 
and so we do not repeat it here. Instead, we 
focus on the illustration of phenomenological 
capabilities of the program that are relevant 
for the description of 6 — u transitions. 

Numerical results reported below are ob- 
tained within the standard framework for per- 
turbative QCD computations. We employ the 
on-shell renormalization for the 5-quark field 
and the 6-quark mass. The strong coupling con- 
stant is renormalized in the MS-scheme. We 
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Figure 1 : The coefficient of the second order correction 
to the lepton invariant mass distribution. The sohd 
curve is the analytic result of Ref . [2^1 . 

note that we do not include the charm mass 
dependence when we compute contributions of 
additional qq pairs to the decay rate. As was 
explicitly shown in Ref. [30| . contributions of 
virtual charm loops for physical value of rric 
can be obtained, with a good accuracy, from 
the bottom quark loops by equating charm and 
bottom masses. We will use this recipe in what 
follows. We write the differential decay rate 
for 6 X^ev through NNLO in perturbative 
QCD as 

dr = dr("' + a.dr'i' + a^ar^^) + o{al), (i) 

where = cks/tt and Us is the MS strong cou- 
pling constant at the scale /i = m?,. By integrat- 
ing the fully differential decay rate over all the 
available phase-space for final state particles, 
we obtain a prediction for the 0{a1) correction 
to the total decay rate. We write the result of 
our numerical integration in the following way 

r(2) = r(°) (-29.98(8) + 2.143(7)iV/ 
-0.0243iV,i) , 

where Nf = 3 denotes the number of mass- 
less quarks in the theory and Nh ~ 2 de- 
notes the number of quarks whose mass co- 
incides with the &-quark mass. Also, F^"'' = 
G%\Vub\^ml/{192n'^) is the total decay rate for 
b — ueP at leading order in perturbative QCD. 
Comparing our computation to the analytic re- 
sults presented in Ref. [1^, we find agreement 
for each term shown in Eq.® to better than 
five per mille. 

Having reproduced the known result for the 
NNLO QCD corrections to the total rate, we 
can now proceed to the discussion of kinematic 
distributions. Our numerical program is set up 
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Figure 2: The cumulative histogram that shows Lqq 
as a function of the cut on the charged lepton energy 
El > Ecut- No cut on the hadronic invariant mass is 
applied. 

in such a way that it can compute various kine- 
matic distributions, both conventional and cu- 
mulative, in a single run. To illustrate this, we 
show in Fig.[T]dr*^^-' /dg^, where q^ is the invari- 
ant mass of the lepton pair. The solid curve 
is the result of the analytic calculation from 
Ref. 28] . The numerical and analytical results 
perfectly agree for all values of q"^ except in the 
region q^ where some discrepancy is ob- 
served. This discrepancy is not surprising since 
the analytic results of Ref. (2^ were obtained 
as an expansion around q^ = ml so that devia- 
tions at small q^ reflect convergence problems 
of the analytic computation in that region. 

To further discuss kinematic distributions, 
we follow Ref. and define moments of the 
partonic invariant mass = (pf, — Pe — PuY 
and energy Ex = Eb ~ ~ E^, in dependence 
of the lower cut on the electron energy i^cut and 
the upper cut on the partonic invariant mass 
Afcut- More specifically, we write 

L,, = {M^x'E'x0iEe~E,,,t)e{M,,,t-Mx)) (3) 

where (...) denotes the normalized phase-space 
average for final-state particles in 6 — >■ XuCD 

We note 20] that one can use Ly's defined in 
Eq. ([3]) to obtain moments of the lepton invari- 
ant mass q^ . 

We write the moments in Eq.Q as an ex- 
pansion in the strong coupling constant and ex- 
plicitly separate the BLM corrections 

(5) 
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Figure 3: The cumulative histogram that shows L^^' 
as a function of the cut on the charged lepton energy 
El > Ecut- No cut on the hadronic invariant mass is 
applied. 



The BLM correction to these moments is ob- 
tained by computing the contributions of a 
massless qq pair i'^^'"^ and then by rescahng 
it by the fuh leading order QCD /^-function for 
three massless flavors 

^(2),BLM^_27/2Lg''"^=\ (6) 

For the numerical calculation of the mo- 
ments reported below, we use = 4.6 GeV 
as the value of the 6-quark pole mass. To 
specify partonic cuts, we introduce the physi- 
cal hadronic invariant mass 



Ml 



A"^ + 2mbAEx + mlM]^, (7) 



where A — mB± — rub = 0.6769 GeV, for our 
choice of the 6-quark mass. We impose a cut 
on Mh and translate it into a cut on Mx and 
Ex using Eq.Q. Throughout the paper, we 
use Mh < 2.5 GeV as the cut on the hadronic 
invariant mass. Also, as we already mentioned, 
the BLM corrections are well-known. They 
were discussed previously in the literature (see 
e.g. Ref. [2^) and, for this reason, we focus on 
non-BLM corrections in the remainder of this 
paper. 

Our results for the moments are presented 

in Figs, [m and in Table ffl In Figs. [IE] we show 

cumulative histograms for the non-BLM con- 

(2) 

tributions to two moments Lqj for j — and 
j = 3, with no cut on hadronic invariant mass. 
In both cases, the x-axis shows the applied cut 
on the the charged lepton energy. These figures 
illustrate that our numerical program works as 
a parton level Monte Carlo integrator and that 
it can reliably compute large number of infra- 
red safe observables for the b X^eD decay 



with various cuts in a single run. This should 
be useful for further studies of charmless decays 
of i?-niesons given, in particular, a large num- 
ber of kinematic cuts employed in experimental 
analyses 

To illustrate the dependence of the non- 
BLM corrections on the applied cuts, we show 
the ratio of non-BLM contributions L^^-' to 
NLO ones L^^\ as a function of the electron 
energy cut in lower panes of Figs. I2I3I We note 
that this ratio is renormalization-scale indepen- 
dent. We are interested in this ratio because, if 
it is independent of Ecut , we could have found 
the corrections to the moments without fully- 
differential NNLO computations. However, it 
is apparent from Figs. [33] that this is not pos- 
sible and that non-BLM corrections have a dif- 
ferent functional dependence on £^cut com- 
pared to the NLO ones. In addition, the cut- 
dependence is strongly moment-dependent and 
it is more pronounced for higher- j moments. 

We will now take a closer look at the nu- 
merical values of the computed corrections. To 
facilitate this, we show in Table [T] our results 
for leading, next-to-leading and next-to-next- 
to-leading order partonic moments Lij com- 
puted with the lepton energy cut of 1 GeV and 
the hadronic energy cut AIh < 2.5 GeV. This 
set of cuts was previously studied in Ref. (20| . 

It follows from Table [T] that similar to the 
total rate BLM and non-BLM corrections have 
opposite size, so that the total result for NNLO 
corrections is smaller than the BLM corrections 
taken alone. The non-BLM corrections seem to 
be more important for lower-moments than for 
higher moments. Indeed, the ratio L,-^"* /L.^j^'^'^ 
decreases monotonically by about a factor of 
1.6, from 0.1842 to 0.112, for z = and j chang- 
ing from J = to J = 3. This trend is also 
visible in the absolute magnitude of the correc- 
tions. Taking as{mi,) — 0.24, we find that for 
i = 0, j — 0, the non-BLM corrections increase 
the moment by about three percent while for 
i = 0,j — 3, they become as small as one per- 
cent. 

While these corrections look small com- 
pared to the current 0(10%) uncertainty in 
the \Vub\ determined from inclusive decays, we 



^We also note that our numerical program is rather 
fast. For example, all numerical results reported in 
this paper, including distributions shown in Figs. 0X2131 
and in Table [T] were obtained in an overnight run on a 
modest-size computer cluster. 
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r(0) 




r{2,BLM) 


r(2) 








0.87135 


-2.261(4) 


-27.7(1) 


5.1(1) 





1 


0.29306 


-0.738(2) 


-8.13(1) 


1.38(2) 





2 


0.10789 


-0.2558(8) 


-2.55(1) 


0.38(1) 





3 


0.04210 


-0.0920(4) 


-0.815(2) 


0.091(6) 


1 





0.0 


0.13110(7) 


2.231(1) 


-0.638(3) 


1 


1 


0.0 


0.05265(3) 


0.882(1) 


-0.256(1) 


1 


2 


0.0 


0.02207(2) 


0.365(1) 


-0.106(1) 


2 





0.0 


4.973(4) -10-^ 


6.83(1) -10-^ 


-9.8(1) -10-^ 


2 


1 


0.0 


2.144(2) -10-^ 


2.93(1) -10-^ 


-4.3(1) -10-^ 


3 





0.0 


3.452(7) •10-'^ 


4.41(1) -10"^ 


-4.9(1) -10-* 



Table 1: Moments of the partonic invariant mass and the partonic energy Ex with the hadronic invariant 

mass cut Mjj < 2.5 GeV and the charged lepton energy cut Ei < 1 GeV. See text for details. 



note that Ref. fi^] estimates the total theoret- 
ical uncertainty on \Vub\ that can be achieved 
with various kinematic cuts on Mh, E^, and 
to be close to six percent. The uncertainty in 
perturbative corrections, which mainly refers 
to non-BLM 0{a^) effects that we discuss in 
this paper, is believed [42*1 to be responsible for 
30 to 50% of the full theory uncertainty. Our 
calculation allows to remove this part of the the- 
ory uncertainty by providing explicit results for 
non-BLM corrections. 

For exam ple, one of the scenarios consid- 
ered in Ref. [42] is a high-cut on the lepton 
energy Ei > 2 GeV; it corresponds to the 
measurement by the BABAR collaboration re- 
ported in Ref. [ill. We find the non-BLM cor- 
rection to Loq{Ei > 2 GeV) using the cumu- 
lative histogram in Fig. [5] and observe that it 
changes the leading order moment L^^^{Ei > 
2 GeV) = 0.257 by 6%I1 Since the experimen- 
tal measurement corresponds to |Kib|^ioOj a 
6% shift in ipo due to non-BLM corrections 
translates into a —3% shift in Vub- We stress 
that the above number is given to illustrate 
the magnitude of the expected effect; a precise 
statement about the impact of non-BLM cor- 
rections requires a dedicated analysis along the 
lines of Ref. 42|. However, it is clear that our 
computation should help in removing a signifi- 
cant fraction of the full theory error in | Vub \ as 
estimated in [i^ for the Ei > 2 GeV cut. 

We also note that it is customary to con- 
sider normalized moments, which are defined 



as Ci 



Lij / Lq 



Since both the numera- 



^For comparison, we note that the corresponding 



BLM corrections to L 



(0) 



-30%. 



tor and the denominator in the definition of 
C'ij receive perturbative corrections, we need 
to consistently expand Cij in a series in as to 
establish how stable it is against radiative cor- 
rections. We find that in case of Cqj, the non- 
BLM corrections are close to one-fifth of the 
BLM corrections for all values of j and they 
change the normalized moment by —0.5% for 
j ~ 1 and by —1.45% for j = 3. 

The situation changes dramatically for par- 
tonic invariant mass moments Lij, with i ^ 0. 
In this case the leading order partonic moments 
vanish since in the b — > ue~P process the par- 
tonic invariant mass is zero. As the result, for 
these moments our NNLO calculation is, es- 
sentially, next-to-leading order and the signif- 
icance of non-BLM corrections increases. As 
follows from Table [T] for Lij moments with i ^ 
and j 0, the non-BLM corrections can be as 
large as 30%. 

To conclude, we presented a computation 
of 0(0^) corrections to the fully-differential 
decay rate of charmless semileptonic b decay, 
b — J> XuSv. Our calculation provides a NNLO 
QCD description of arbitrary infra-red safe ob- 
servables and allows arbitrary kinematic cuts 
including those that closely match the ones 
employed in experimental analyses. We con- 
structed a parton-level Monte-Carlo integrator 
which can be used to compute large number 
of relevant observables and kinematic distribu- 
tions in a single run of the program. This calcu- 
lation, together with earlier results on NNLO 
QCD corrections to fully-differential b — > civ 
transition 44, 4^, makes all inclusive semilep- 
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tonic decays of 6-quarks upgraded to that accu- 
racy. We hope that these resuhs wih contribute 
to the reduction of the theoretical error on \ Vub\ 
and \ Vcb\ that will be achieved in the forthcom- 
ing _B-physics experiments. 
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